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Abstract 

The question of whether layer decouphng and vortex-lattice melting occur separately or 
not in the mixed phase of pristine layered superconductors in the extreme type-II limit is 
studied through a partial duality analysis of the layered XY model with uniform frustra- 
tion. We find that both transitions occur simultaneously if the normal/superconducting 
transition of the vortex lattice in an isolated layer is first order and if a sufficient degree of 
layer anisotropy exists. We also find that a crossover to a highly entangled vortex lattice 
phase with relatively low phase rigidity across layers does not occur in practice under any 
circumstances at temperatures below the two-dimensional vortex-lattice melting transition. 

PACS Indices: 74.60.Ge, 74.60.-w, 74.60.Ec, 74.25. Dw 
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I. Introduction 



It is now well established experimentally that the Abrikosov vortex lattice state in 
clean high-temperature superconductors undergoes a first-order melting transition into a 
liquid phase. High-temperature superconductors are layered and extremely type-II.^ The 
former vortex liquid phase in the most anisotropic materials like BSCCO is best described 
by a liquid of planar vortices inside of decoupled layers.^ A longstanding question is whether 
melting and layer decoupling occur simultaneously as a sublimation transition, or whether 
a separate decoupling transition follows the melting transition. Some experimental studies 
on the highly anisotropic BSCCO material show evidence for sublimation,^'^ while most 
experimental studies of the less anisotropic YBCO material point to separate melting and 
decoupling transitions.^ 

The experimental situation outlined above suggests that the degree of anisotropy is 
what in fact determines whether or not the vortex lattice in a layered superconductor 
sublimates. We shall study this proposal theoretically by analyzing the layered XY model 
with uniform frustration, which provides a qualitatively correct description of the thermo- 
dynamics deep inside of the mixed phase in extremely type-II layered superconductors.^ 
After performing a partial duality transformation on the XY model that is particularly 
well suited to the weak-coupling limit, ^ we find that there can exist as many as three 
different decoupling transitions at temperatures Td < Tm < Tx , respectively. (We use the 
term 'transition' here loosely to describe both genuine phase transitions and cross-overs.) 
The phase correlation length across layers is equal to the inter-layer spacing along the 
dimensional crossover line^'^° at T = Tx that separates two-dimensional (2D) from three- 
dimensional (3D) vortex- liquid behavior.^ The phase correlation length across layers then 
either diverges or jumps to infinity along the melting line, T = Tm, which separates the 
superconducting and normal phases. Last, the crossover line T = Td that lies inside of 
the ordered phase is defined by the point at which the Joscphson coupling energy reaches 
about half of its zero-temperature value. The macroscopic phase rigidity across layers be- 
comes small in comparison to its zero-temperature value at temperatures T > Td because 
of the entanglement of fluxlines between adjacent layer s.®'-^-*^ All three decoupling transi- 
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tions occur separately in the continuum regime at low perpendicular vortex density, but 
Td crosses below the 2D melting temperature at only exponentially weak inter-layer cou- 
pling. At a moderate concentration of vortices, on the other hand, we find that the three 
decoupling transitions collapse onto a single sublimation line for weak enough Josephson 
coupling. This is due to the first-order nature of the melting of the 2D vortex lattice 
in such case. These results are compared with previous theoretical calculations based on 
the elastic medium description of the vortex lattice^'^'-^-*^'-^^ and with direct Monte Carlo 
simulation results of the XY model itself.^ 

II. Duality Theory 

The layered XY model with uniform frustration is the minimum theoretical descrip- 
tion of vortex matter in extremely type-II layered superconductors. Both fiuctuations of the 
magnetic induction and of the magnitude of the superconducting order parameter are ne- 
glected within this approximation. The model hence is valid deep inside the interior of the 
mixed phase. The thermodynamics of the 3D XY model with anisotropy and uniform frus- 
tration is determined by the superfluid kinetic energy E^y = ~ J2r ^. Jij.^^^[^ij.4> ~ ^/x]|r, 
which is a functional of the superconducting phase (f){r) over the cubic lattice. Here, 
Jx — J — Jy and Jz = are the local phase rigidities, with anisotropy parame- 

ter 7' > 1. The vector potential = (0, 27r/a;/a, 0) represents the magnetic induction 
oriented perpendicular to the layers, B± = ^of/a^. Here a denotes the square lattice 
constant, which is of order the zero-temperature coherence length, $0 denotes the flirx 
quantum, and / denotes the concentration of vortices per site. The component of the 
magnetic induction parallel to the layers in taken to be null throughout. 

We shall now analyze the above layered system in the selective high-temperature 
limit, ^ J^. Following ref. 8, the corresponding high-temperature expansion can be 
achieved through a partial duality transformation of the layered XY model along the z axis 
perpendicular to the layers. This leads to a useful layered Coulomb gas (CG) ensemble in 
terms of loops of Josephson vortices in between layers (fiuxons).^"^ In particular, suppose 
that I denotes the layer index, that r represents the x-y coordinates, and that r = (f, I). 
Phase correlations across N layers are then described by the phase auto-correlation function 
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probed at sites set by an integer field p{r) = 5r",o((^/,i — (^/,iv)- These can be computed from 
the quotient 

(exp[i5^p(r)(/>(r)]) = ZcgM/^cg[0] (1) 

r 

of partition functions for a layered CG ensemble that describes the nature of the Josephson 
coupling:^ 

ZcgIp]^ J2 yo^'''^-^iC[qi]-e-'^r-^^% (2) 

{n.(r)} 

where nz{f, /) is an integer field on links between adjacent layers I and I + 1 located at 2D 
points r. The ensemble is weighted by a product of phase auto-correlation functions 

C[qi] = (exp [i J2 Qiimr, O] )j (3) 
f 

for isolated layers I probed at the dual charge that accumulates onto that layer: 

qiir) = p{f, I) + n,{r, I - 1) - n,{r, I). (4) 

It is also weighted by a bare fugacity yo that is raised to the power N[nz\ equal to the 
total number of dual charges, rig = ±1. The fugacity is given by yo = Jz/'^ksT in the 
selective high-temperature regime, <^ ksT, reached at large model anisotropy. Also, the 
average number of riz charges per link is equal to^ 2yo((cos — yo), which is less than 
Jz/ksT. This implies that the layered CG ensemble (2) is dilute in such case, because 
j/o ^ 1- The former is required by the approximate nature of Eq. (2), which neglects 
multiple occuppancy of the dual charges, n^, on a given link. Last, the thermodynamics of 
the layered XY model is encoded by its partition function, which is given by the following 
product: 

4'y[0] = [Wz/kBT)]^' ■ Zcg[0] • niZ^iln. (5) 

Here, Io{x) is a modified Bessel function, and [0] partition function of an isolated 

layer. Also, M' denotes the total number of links between adjacent layers. 

Interlayer correlations of the layered XY are easily determined using the CG ensemble 
(2) when the phase correlations within an isolated layer are short range. ^ Let us introduce 
the notation (f)i,i'{r) = (f){f,l') — (j>{r,l) and take = due to the null magnetic field 



parallel to the layers. A useful (in)equality for the autocorrelator between any number of 
layers, n + 1, can be computed to lowest order in the fugacity, j/o- It reads* 



where 



Co 



(27r)2\Co 



{yoCo/a'r, (6) 



C,= / d^r|C(r)|e*^- (7) 



is the Fourier transform of the magnitude of the phase auto-correlation function (3) for an 
isolated layer probed at two points, fi and f2'- 

C{l,2) = \C{r,2)\e-'I^^'^'^-'', (8) 

where A' is a suitably gauge-transformed vector potential (see below). Its magnitude de- 
pends only on the separation fi2 = ri—f2 between the probes, and it decays exponentially 
at separations beyond a characteristic correlation length $^2D due to the phase-incoherent 
state that is presently assumed. The layered CG ensemble (2) is therefore in a confining 
phase. The prefactor in brackets above in Eq. (6) typically decays polynomially with 
the separation n between layers. Also, Eq. (6) is an equality for n = 1,^^ as well as for 
pure gauges such that A' = V0o (see below). To conclude, the autocorrelator (e*'^'''+") 
across layers decays at least exponentially with the separation n in the weak-coupling limit, 
yo 0, of the disordered phase. 

The layered CG ensemble (2) can also be used to determine interlayer correlations 
in the ordered phase. Consider again an isolated layer, and suppose that general phase 
auto-correlation functions (3) are quasi-long range: 

CM = 9o^ ■ exp V2D ?(l)ln(ri2/ro) ?(2)] • exp(^z ^ 5 • 0o) , (9) 

(1,2) 

where Qq is equal to the phase rigidity of an isolated layer in units of J, where ?i+ is equal to 
half the number probes, where tq is the natural ultraviolet scale of order the inter-vortex 
spacing, = a/f^^'^, and where (j)o{r) should resemble the unique zero-temperature 
configuration (independent of the layer index, I). The system of dual (rig) charges in the 



layered CG ensemble (2) is then in a plasma phase at low temperatures r]2D < 2.^'-^^ In 
such case, the macroscopic phase rigidity across layers is approximately given by^ 

p^/J,^(cos0^+i)-yo. (10) 

Furtheremore, in this case an appropriate Hubbard-Stratonovich transformation of the 
CG partition function (2) in the absence of a source {p = 0) reveals that it is equal to the 
corresponding one Zld[0] = / T>0 e~^^°^'^^'^ for a renormalized Lawrence-Doniach (LD) 
model up to a factor that is independent of the Josephson coupling, J^. The corresponding 
energy functional is given by^ 
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where J = kBT/27rr]2D is the macroscopic phase rigidity of an isolated layer, ^® and where 
Ao = Ya is the Josephson penetration length. The above continuum description (11) is 
understood to have an ultraviolet cut off of order the inter-vortex spacing, tq. A standard 
analysis of the product of partition functions (5) then yields that the strength of the local 
Josephson coupling is given by 

(cos0i,/+i) = yo +5'o(cos6'«,/+i), (12) 

where 9i^i+i = 9i+i - 9i. 

To compute (cos^^^z+i) in the weak-coupling limit, it is sufficient to consider only 
layers I and Z + 1 in isolation from the rest of the system. At low temperature r]2D <S 1, the 
harmonic approximation for the Josephson coupling term in Eq. (11) is valid: cos9i^i+i = 
1 — ^9fl_^-^^. The resulting gaussian integration then yields (cos = e~^^'''+i^/^, with 
i^il+i) ~ ln(Aj/ro). Here Aj is of order the Josephson pentration length, Aq = 7'a. 
Subsitution into Eq. (12) then produces the result^ 

(cos0^+i) = yo + 5^o(ro/Aj)^=^ (13) 



for the strength of the local Josephson coupling at low temperature r]2D ^ 1- The latter 
agrees with the result produced by analyzing a fermion analogy for the LD model (11), 
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as well as with an estimate by Glazman and Koshelev for the zero-field case (ro ~ a). 
Subsitution of this result into Eq. (10) therefore yields the formula 



for the macroscopic phase rigidity across layers in this regime. To conclude, macroscopic 
phase coherence exists across layers in the ordered phase (9). 



We shall now review the phase diagram that results from employing the above duality 
analysis for the layered XY model in the continuum limit, ^ a — 0, which coincides with the 
regime of small perpendicular flux density, / <^ 1/36. In the absence of surface barriers, 
Monte Carlo simulations^^ indicate that the vortex liquid phase of an isolated layer solidifies 
into a "fioating" vortex lattice phase at the 2D melting temperature, ksTm^^ = J/20. A 
recent duality analysis of such a single layer finds that the standard 2D melting scenario-*^^ 
takes place as long as rigid translations of the 2D vortex lattice are prohibited by surface 
barriers. In particular, general phase auto-correlation functions follow the form (9) in 
the vortex lattice phase at T < T^^\ with a 2D correlation exponent that takes on an 
extremely small value-*^^ 'q2D — (287r)~"'^ just below the 2D melting temperature, Tm^\ 
Further, r}2D decreases linearly to zero with decreasing temperature in the 2D vortex 
lattice. On the otherhand, the phase auto-correlation function (8) decays exponentially 
with separation in the hexatic phase that lies at temperatures just above Tm^\ The 
associated correlation length, ^2D, diverges exponentially as temperature cools down to 
Tm^\ The auto-correlation function retains, however, the trivial phase factor of the 2D 
vortex lattice:^^ jl A' ■ dr = 0o(2) - 0o(l). 

We now illustrate that there exist as many as three distinct decoupling temperatures:^ 
Tx > Tm > Td- Consider the weak-coupling limit of the layered XY model, 7' 00. Eq. 
(6) then becomes an equality in the hexatic phase of an isolated layer due to the trivial 
phase factor in the phase auto-correlation function (S).-*^^ The phase correlation length 
across layers, is therefore equal to the spacing d between adjacent layers when 



Ps/Jz = 9o{ro/Aj) 



(14) 



III. Continuum Limit 




(15) 
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This defines a dimensional cross-over field,^'^ ^° 

/7x ~ 9o{J/kBT){i2D/a^^f (16) 

in units of the naive decouphng scale ^q/Aq, that separates 2D from 3D vortex- liquid 
behavior.^ It is traced out in Fig. 1. In these units, to be used hereafter, /7'^ gives the 
perpendicular field. The system is best described by a decoupled stack of 2D vortex liquids 
at fields above /7x . On the ordered side at T < Tm^\ Eq. (14) for pj- implies that long- 
range order across layers exists: ^± = oo. And since go J is equal to the phase rigidity 
of an isolated layer, Eq. (14) also implies that 3D scaling is violated at weak-coupling, 
(ro/Aj)*'2o <^ 1^ in which case the phase rigidity across layers, p^, is small in comparison 
to its value at zero temperature, Jz- This occurs at fields above the decoupling scale 

= however, which is astronomically large and of order 10"^^ at temperature 

below 2D melting due to the extremely small bound on the correlation exponent there, 
r]2D < (287r)~-'^. At large anisotropy, 7' > 7^) the system is best described by an entangled 
stack of 2D vortex lattices-*^-*^ that exhibit a relatively small macroscopic Josephson effect.^ 
Last, the CG ensemble (2) indicates that a 3D vortex-lattice melting transition occurs 
at an intermediate temperature Tm when the typical distance between neighboring dual 
charges, Ug — ±1, grows to be of order ^20, at which point these charges are confined 
into neutral pairs. ^'^^ It can be shown that Tm lies inside of the 2D-3D cross-over window 
[Tm^\Tx] by virtue of this definition (see ref. 8, Eq. 62). Also, by comparison with 
the layered CG ensemble (2) in zero field, ^'^^ the author has argued that in the weak- 
coupling limit, Tm marks the location of a second-order melting transition that separates 
the superconducting and normal phases. This means that ^_l(T) diverges as T cools down 
to Tm- A second-order transition in the vortex- liquid phase of YBCO that resembles the 
above has been reported recently.^ 

Let us now determine what happens as interlayer coupling increases from the weak- 
coupling limit just studied. The charges are screened at low temperature, T < Tm^\ 
which means that no phase transition can take place as a function of the anisotropy 
parameter,^ 7'. Instead, a cross-over region exists for anisotropy parameters below 7^ 
that separates a set of weakly coupled 2D vortex lattices at high field from a conventional 
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3D vortex lattice at low field. Again, the extremely small bound on the 2D correlation 
exponent 7720 at temperatures below 2D melting indicates that the former weakly coupled 
phase is not attainable there in practice. Eqs. (13) and (14) also imply that the Joseph- 
son effect is essentially independent of field/anisotropy at these temperatures, T < Tm^\ 
This observation is consistent with Monte Carlo simulation results of the layered XY model 
with uniform frustration.^ On the disordered side, T > Tm^\ the phase correlation length 
across layers, , begins to grow larger than the spacing between adjacent layers at fields 
below /7x^.^ Outside of the 2D critical region, at ^20 ~ ^vx? Monte Carlo simulations of the 
layered XY model with uniform frustration indicate that first-order melting occurs along 
the decoupling contour (cos0i,/+i) ~ 1/2.^'^^ The resulting phase diagram is depicted by 
Fig. 1. 

IV. Sublimated Decoupling 

We shall next apply the partial duality analysis outlined in section II to the layered XY 
model with only moderately small frustration. Let us consider again an isolated XY model 
over the square lattice, but with a uniform vorticity (frustration) between l/30</<l/2. 
Monte Carlo simulations indicate that a depinning transition at ksT^^^ = 1.5/J now 
separates a pinned triangular vortex lattice at low-temperature from a vortex liquid phase 
at high temperature.^^ The depinning transition is first order and no signs of a "fioating" 
vortex-lattice phase are observed. Strict long-range phase correlations then exist at low 
temperatures T < T^^^ in the pinned phase, which implies that the phase auto-correlation 
functions are given asymptotically by Eq. (9) with ri2D = 0. Also, the disordered phase 
at high temperature T > Tp"^^^ should be hexatic due to the underlying square-lattice 
grid.^^ This means that the phase autocorrelations (8) exhibit exponential decay as well 
as a trivial phase factor: ^20 < 00 and • = 0o(2) — 0o(l)- Below, we shall use 

these facts to map out the phase diagram of the layered XY model at such relatively high 
vortex density. 

The first-order nature of the depinning transition in an isolated XY layer with rela- 
tively large uniform vorticity, 1/2 > / > 1/30, implies that the phase correlation length 
is finite at temperatures just above the depinning transition: ^2DiTp'^^^+) < 00. By Eq. 
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(16), the 2D-3D cross-over field here must also then be finite. Notice that f^'^ is larger 
than unity at depinning if > and if ^ro ~ 1; since J > ksT^^^^ for / < 1/2. Strict 
long-range phase coherence (?72D = 0) exists on the low-temperature side at T < T^^\ 
however. We therefore reach the remarkable conclusion that at large anisotropy parame- 
ters of the corresponding layered XY model, 7' ^ Txl^p^^"*]; ^i^® ^ — Tp^^^ marks 
a sublimation transition that separates a decoupled vortex liquid at T > T^^^ with es- 
sentially no interlayer phase coherence, ^± < ci, from a pinned 3D vortex lattice state at 
T < Tp'^^^ with long-range interlayer phase coherence, = 00. As depicted by Fig. 2, 
no 2D-3D cross-over regime exists in such case. Also, comparison of Eqs. (13) and (14) 
with the fact that the 2D correlation exponent 7720 vanishes in the low-temperature phase 
implies that the cross-over at 7' = 7^ (T) between weakly coupled and moderately coupled 
vortex lattices must collapse onto the depinning line at T = Tp"^^^ and 7' > 7x[^p^^^]- 
Indeed, Eq. (13) indicates that the Josephson coupling {cos (pi^i+i) is independent of field, 
/7'^, at temperatures below the sublimation transition and at such large anisotropy pa- 
rameters. Last, the local Josephson coupling jumps down to a small value given by the 
vortex- liquid result, -"^^ Eq. (6) at n = 1, once the vortex lattice sublimates. Similar jumps 
of order unity have been observed at vortex-lattice melting in BSCCO.^ In conclusion, the 
three possible decoupling transitions collapse onto a single sublimation transition! Such 
point-like as opposed to line-like melting of the vortex lattice has been observed in Monte 
Carlo simulations of the layered XY model with moderately small frustration.^ 

V. Discussion and Conclusions 

Among the important theoretical results listed above is the local Josephson cou- 
pling in the vortex-lattice phase, Eq. (13), which can be expressed as (cos0/,i+i) = 
yo + 9oe~^'^/^°^^^\ with a temperature scale kBTD{B±) = 2TTj/\n{B^/B\). Here, B\ = 
^»o/Ao is the naive decoupling field^ and J = kBT/27rr]2D is the 2D phase rigidity. -"^^ As ob- 
served previously, the weak logarithmic field dependence above implies that (cos^; is 
of order unity at low temperatures T < T^^^ and at perpendicular fields below the astro- 
nomically large scale Hr, ~ lO^^S^. The local Josephson coupling (13) shows essentially 
no field dependence in such case. This is confirmed directly by Monte Carlo simulations of 
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the layered XY with low uniform frustration/ Despite the fact that the decoupled vortex- 
lattice state characterized by a small "cosine" does not exist in practice at temperatures 
below 2D melting, it is nevertheless remarkable that Td{B±) coincides, to within a large 
numerical constant, with the temperature scale for layer decoupling induced by the un- 
binding of topological defects of the vortex lattice known as "quartets" . These consist of 
two opposing dislocation pairs in parallel inside of a given layer. ^'^^ Comparison with the 
present results then indicates that layer decoupling is indeed due to such a "quartet" un- 
binding mechanism, but that this occurs only for exponentially weak Josephson coupling at 
temperatures below 2D melting (cf. ref. 12). Glazman and Koshelev have also calculated 
the local Josephson coupling (cos0i,/+i) within the 3D elastic medium description for the 
vortex lattice,^ where they find a much stronger dependence T^(5j^) ~ [B*^/ B<^Y/'^Tm^^ 
for the decoupling temperature scale with field, on the other hand.^ This discrepancy 
is due to the fact that the elastic-medium approximation represents a continuum the- 
ory. It therefore accounts only for long-wavelength fluctuations of the phase difference 
across layers. In the weak-coupling limit, the dominant contribution to the "cosine" is 
due to short-wavelength phase fluctuations between adjacent layers. These fluctuations 
are missed by the 3D elastic medium approximation, and we believe that this is why the 
Glazman-Koshelev result^ underestimates the size of the decoupling temperature scale at 
weak coupling. 

In conclusion, a partial duality analysis of the layered XY model with uniform frus- 
tration flnds that sublimated melting/decoupling of the 3D vortex lattice occurs if (z) the 
superconducting- normal transition of an isolated layer is first-order and if (m) a sufficient 
degree of layer anisotropy exists. Condition (i) is gauranteed at strong substrate pinning, -"^^ 
1/2 >/> 1/30. It has also been emphasized that no decoupled vortex-lattice state ex- 
ists at temperatures below 2D ordering except for exponentially weak Josephson coupling 
between layers (see Figs. 1 and 2). This is notably consistent with complementary calcula- 
tions that include interlayer magnetic coupling, but that turn off the Josephson coupling. 
It must be mentioned, however, that the magnetic coupling between layers is weak in the 
extreme type-II regime studied here, and that this coupling can in fact be incorporated 
into the present duality analysis (2) of the vortex lattice in layered superconductors via an 
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effective "substrate potential" for isolated layers (see ref. 12). The additional substrate 
consists of an array of commensurate pins that mimics the magnetic effect of the vortex 
lattice in adjacent layers. It can therefore only increase phase coherence (3) inside of each 
2D vortex lattice. This means that the bound, rj2D < (287r)~^, on the phase correlation 
exponent of the 2D vortex lattice continues to hold. Hence, within the "substrate po- 
tential" approximation for magnetic coupling, -"^^ the decoupling crossover to an entangled 
vortex lattice^^ with <^ does not occur in practice at temperatures below 2D melting 
in the extreme type-II regime [see Eq. (14)]. We remind the reader that rigid translations 
of the vortex lattice are assumed throughout to be prohibited by surface barriers (see ref. 
19). 

The author is grateful for the hospitality of the Instituto de Ciencias de Materiales 
de Madrid, where this work was completed, and to Marty Maley and Paco Guinea for 
discussions. 
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Figure Captions 



Fig. 1. Shown is the proposed phase diagram for the layered XY model with uniform frus- 
tration in the continuum regime, / ^ 1/36. Notice the absence (in practice) of a 
decoupling transition at temperatures below 2D melting. Rigid translations of the 
vortex lattice are assumed to be prohibited by surface barriers. The mean- field tem- 
perature dependence J oc T^o — T is also assumed. 

Fig. 2. The proposed phase diagram for the layered XY model with moderate uniform frus- 
tration, 1/30 < / < 1/2 is displayed. The mean- field temperature dependence 
J (X Tco — T is assumed once again. 
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Erratum: "Sublimated decoupling of the vortex lattice 
in extremely type-II layered superconductors", 
[Phys. Rev. B 66, 214506 (2002)] 

J. P. Rodriguez 

The decoupling field for temperatures that lie below the 2D ordering transition that 
was derived in the discussion following Eq. (16) is more generally given by 

/7g = (ro/avx)'eV^-, 

where tq ~ ayx was implicitly assumed. Although the latter is not necessarily true, the 
ratio ro/ovx must be larger than = We have k ~ 100 in YBCO for example. 

The above then implies that the decoupling field is bounded by > 10^^ at temperatures 
below 2D ordering in such case, since ri2D < (287r)~^. It therefore remains exponentially 
big. 

More seriously, the claim made in section IV that the 2D phase correlation exponent is 
null at temperatures that lie below the 2D vortex-lattice depinning transition is incorrect. 
What is null is its vortex component, which leaves the spin-wave result 7)20 = kBT/2'KJ 
for the net exponent. The sentences in the middle of both paragraphs of section IV that 
begin with "Strict long-range phase ..." must therefore be replaced with "Quasi long-range 
phase Also, the equation "7720 = 0" that appears in both of these sentences must be 
replaced with "7720 = ksT /2ttT\ The rest of section IV remains valid for Josephson 
coupling that is not exponentially weak. The equation displayed above, for example, 
yields an astronomically large lower bound > (ro/avx)^ • 10^^ for the decoupling field 
at temperatures below 2D ordering and at an in-plane vortex concentration of / = 1/25. 
This bound is due to the value ksTp'^^ = 0.06J of the first-order transition temperature 
of an isolated layer in such case. 

The above corrections do not change any of the conclusions drawn in the paper. 
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